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Abstract 

We study a simplified model of the electronic structure of compounds of the 
type of Lai_a;Sr x Mn03. The model represents each Mn 4+ ion by a spin S=l/2, 
on which an electron can be added to produce Mn 3+ . We include two strong 
intratomic interactions in the Hamiltonian: exchange (J) and Coulomb (£/). 
Finally, to represent the effect of Sr substitution by La in a simple way, we 
include a distribution of diagonal energies at the Mn sites. Then we use Green 
function techniques to calculate a mobility edge and the average density of 
states. We find that according to the amount of disorder and to the concen- 
tration of electrons in the system, the Fermi level can cross the mobility edge 
to produce a metal to insulator transition as the magnetization decreases (in- 
crease of temperature). If the disorder is large, the system remains insulating 
for all concentrations. Concentrations near zero or one favor the insulating 
state while intermediate values of concentration favor the metallic state. 
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The recent discovery of extremely large magneto resistance in Lai_. r A. r Mn03 (where A 
stands for Ca, Sr, and Ba) and other related oxides U has renewed interest in these types 
of compounds. The magneto resistance values found at large fields seem to be connected 
with a metal-insulator transition found in some of the compounds at temperatures equal or 
lower than the ferromagnetic to paramagnetic transition. The T vs. x phase diagram of 
these compounds is quite similar to that of the high temperature superconductors at low 
and intermediate concentrations of A, but with the ferromagnetic phase substituting the 
superconducting one. At the end of the composition range both LaMnO"3 and CaMnO"3 are 
antiferromagnetic insulators, while in the intermediate composition range the compounds are 
ferromagnetic metals. Much has been learned about these compounds and alloys since the 
earlier studies of Jonker and Van Santen ||. The later experimental results including single 
crystal measurements are listed in Ref. 3. Zener [|J advanced a mechanism to understand 
the properties of the different compounds. Based on the coexistence of Mn 3+ (3d 4 ) and Mn 4+ 
(3d 3 ) within the doped materials, he proposed that the displacement of the fourth electron 
between them produces the metallic conductivity and at the same time provides a mechanism 
for the ferromagnetic alignment of the spins. Anderson and Hasegawa || calculated the 
interaction between two magnetic ions mediated by the " double exchange mechanism" . They 
point out the importance of the orbital degeneracy for this mechanism as is the case in the 
Mn ions, de Gennes || made a study of the competition between superexchange and double 
exchange in the dilute limit (x « 1) compounds and proposed the existence of canted or 
spiral magnetic structures. Kubo and Ohata used a spin wave approach to study the 
temperature dependence of the resistivity well below the Curie Temperature and a mean field 
approach to the many body Hamiltonian to obtain the electronic and magnetic excitations 
in the metallic phase. Mazzaferro, Balseiro, and Alascio || used a strong coupling approach 
based on the similarity of these compounds to the intermediate valence Tm compounds. 
They introduced in the model the possibility of non- equivalent sites for the Mn ions. To 
our knowledge, this was the first time that the possibility of a metal to insulator transition 
was mentioned in connection with double exchange. Since magnetization and conductivity 



are intimately connected it is evident that any appropriate description of the properties of 
these compounds must include those factors that affect the mobility of the carriers. 

Furukawa used the infinite dimension Kondo lattice model with classical spins to describe 
several properties of the compounds. Within the model, he was able to fit the resistance 
versus magnetization curves of Tokura et. al. showing that double exchange is essential to 
the theory of doped LaMnOs. Millis et. al. |10| analyze a 'ferromagnetic Kondo Hamiltonian' 
using different approximate methods to conclude that double exchange alone is not sufficient 
to describe the properties of Lai-^Sr^MnOs. 

In this paper we extend the study of Ref. 8 including a continuous distribution of 
inequivalent sites for the Mn ions. We find that as the magnetization increases from zero 
to saturation the system may undergo a transition from a state where the Fermi level falls 
below the localization edge to one where it falls above the edge thus substantially changing 
the transport properties. 

We consider a model were one electron can hop from the 3d 4 configuration of Mn 3+ to 
the 3d 3 configuration of a Mn 4+ nearest neighbor. To consider the Zener double exchange 
mechanism, we include an exchange energy between the fourth electron and a core formed 
by the three localized electrons in the lower d orbitals at each site. To avoid double occu- 
pation by itinerant electrons, we take a large Coulomb interaction between them. Finally 
to simulate the effect of disorder produced by the substitution of La by Sr or other divalent 
ions and other defects, we introduce a distribution of diagonal energies for different sites. 
We use a very simplified model Hamiltonian that represents each Mn 4+ ion at site i by a 
spin Si, on which one electron can be added to produce Mn 3+ (for simplicity in what follows 
we take Si = |). When an electron is added in the d-shell of site i, an exchange coupling 
J is included to favor parallel alignment of the added electron to the already existing spin 
[ ITJ . Also to avoid the possibility of Mn 2+ we include a strong Coulomb repulsion U. The 
Hamiltonian reads: 



i,a <i,j>,<7 



where cj CT , c ia creates and destroys an itinerant electron with spin a at site i, respectively. 
Si and <Tj are the Pauli matrices for spin | at site i for localized and itinerant electrons 
respectively, q is the on-site energy, t the hopping parameter between nearest neighbors, U 
the on-site Coulomb repulsion between two itinerant electrons, and J is the ferromagnetic 
(J > 0) coupling between the localized and itinerant electrons. As is well known, this 
coupling does not give rise to divergence in the impurity scattering as in the Kondo case nor 
to heavy electron dynamics as in the Kondo compounds. Without losing essential physics 
we simplify further by taking only the z component of the exchange interaction. Thus the 
states of the system are characterized by itinerant electrons moving on a frozen distribution 
of localized up or down spins. To obtain site Green functions and thus local density of states 
for this problem, we start from Hamiltonian (1) for t = and we calculate the local Green 
function. For spin up we obtain: 

G^ a -« c lh c A »- {u) _ Eia){u _ Eia _ u) ■ (2) 

where =< cj^c^ >, a = + (a = — ) for up (down) localized spin, and E ia = o L — aJ. In 

order to obtain an approximate solution to the problem for (t ^ 0) at the start we ignore 

the site dependence of the diagonal energies: i.e. we set e% = e and we are left with a binary 

alloy problem. Using the Renormalized Perturbation Expansion (RPE) ]I2] and Eq. (2) for 

U — ► oo we obtain the corresponding local Green functions. Spin up gives 

_ (1-ng) 
ira [co - E ia - (1 - %)A lT ] ' [6) 

where Aj-j is the corresponding self-energy given by 

A,T = t 2 E G' l+ s 1a , (4) 

6a 

the summation over 5 runs over the nearest neighbors to site i and the G' are the local 
propagators avoiding i. Their self-energies A- +(5 | are given by 



A'i+si —t 2 ^2 G" +S+S ,^ a , (5) 

S'a 

and so on. Thus, for a system involving an infinite number of sites the RPE iterated proce- 
dure as indicated above gives an infinite number of steps. In order to get an approximate 
solution we average over spin configurations the second term in Eqs. (4) and (5) which gives 



Af = (K + l)i 



(1-n+K 



[u- E + -(l-n+)A\ 



(l-nT)v- 



+ [w-£;_-(l-fij-)A' r ]'' (6) 
where E a = (e — aJ), (nj) is the average number of itinerant electrons with spin down 
at localized spin up (down) sites, (K + 1) is the number of nearest neighbors (six for the 
simple cubic Mn lattice), v + (z/_) is the probability that a site has up (down) localized spin 
(notice that (v + + vJ) = 1), and 

(1 — nj)v~ 

+ [u-E_-(l-fij)A\y ' (7) 
Notice that within this approximation the iterative procedure closes in the last Equation. 
This allows us to obtain average density of states for spin up and down. At J > \f~Kt the 
density of states splits into two bands centered at E a with weights and widths that depend 
on the number of sites with each spin. i.e. they depend on the magnetization of the system. 
In Fig. 1 we show the densities of states for spin up and down for several values of the 
magnetization (m), with m = (v + — v_). For J » VKt, the term proportional to z/_ can 
be neglected in Eq.(6) and (7). Similarly for Aj and Aj, consequently it is easy to see that 
— > and — > 0. Eq.(7) reduces to 



A , = 0^) ±i /(^)!_^ +i (8) 



where E = (e — J) and Af results in this case 

(K + l)t 2 ^ 



(9) 



Finally we obtain 

_ (K - - E) T (k + l)y/(u - E) 2 - 4Kt 2 u + 
T+ 2[{K + l) 2 t 2 v + - (u - E) 2 ] ' { ] 

Eq. (8) allow us to obtain the density of states per site as 



(K + l)y4Kt 2 v + — (u — E) 2 
P0T+H = 2n\(K + l)H*v + - (u - E)*\ ■ (H) 

At this point, we introduce the site dependent diagonal energies. As is well known, since 



Anderson's original paper [13] , a distribution of diagonal energies produces localization of 
the electronic states from the edges of the bands to an energy within them which is called 
"mobility edge" (ME). The precise position of the ME is difficult to calculate, different 



localization criteria result in different values for it [F3|. Our concern here is with the changes 
of the ME with magnetization, which do not differ much among the different criteria. For 
this reason we report results using Anderson's original criterium. 

Since the resulting structure of localized and extended states in each band can only 
be qualitatively obtained after having tried several distribution of energies, we present the 
simplest mathematical model for the distributions of energies: i.e. a Lorentzian j0| distri- 



bution of width T. The ensemble-averaged Green function allows us to write an approximate 
density of states given by 

/ + 0O 
PoT+(e')£(e - e')de' , (12) 
-oo 

where L(x) is a Lorentz distribution given by 

L(x) = k/+ n\ • (13) 

In a similar manner we obtain 

In Fig. 2 we show the resulting densities of states and mobility edges for different 
values of the magnetization. Inspection of Fig. 2 shows that the system can change from 
metallic at full magnetization (the Fermi level falls above the mobility edge of the majority 
band) to localized at the paramagnetic state (the Fermi level below the mobility edge of 



the same band). Depending on the values of T and n , situations where the system remains 
localized or metallic for all values of the magnetization can also be found. This possibility 
of the system of changing character of states at the Fermi level from extended to localized 
can be summarized in the phase diagram depicted in Fig.3, where the localized or extended 
character of the states at the Fermi level are shown as functions of the electron concentration 
and magnetization for different values of Y. 

We do not derive here an expression for the Free energy of the system that would allow 
us to calculate different thermodynamic properties, however we can easily estimate the effect 
of disorder on the Curie temperatures: The energy difference between the zero-temperature, 
fully polarized system and the paramagnetic state can be calculated from the density of 
states and is a function of the disorder parameter T and the bandwidth. For T >> \fKt 
the total density of states (up and down bands) is weakly dependent on m while in the 
opposite case the density of states changes substantially as m decreases from saturation 
to zero. The entropy change between paramagnetic and ferromagnetic states, however, is 
practically independent of the parameters in the Hamiltonian and is dominated by the spin 
entropy and approximately equal to ksln^i + n). In fig 4. we show the electronic energy 
as a function of the n for two values of the disorder parameter. The Curie temperatures 
estimated using the above argument range from OK to about 800K at the top of the scale. 
One can see that the Curie temperatures decrease with increasing Y. 

The absence of Mn 2+ and the value of the saturation magnetization indicate that our 
assumptions U, J » t,T are valid in most of the Mn oxides under consideration. The 
low energy properties then, depend only on t, n, and T. Of these, T is the most difficult 
parameter to estimate, as it should include all effects of disorder. Among these, substitution 
of La by divalent ions and O vacancies are the most evident sources of disorder, but other 
types of defects also contribute to T. Polaronic effects derived from Jahn- Teller as well as 
breathing phonons should be present in these materials. It is conceivable that they could 
also be represented within our assumptions of a distribution of diagonal energies. 

We propose here a model that allows to understand the properties of the different magne- 



toresistive compounds in terms of a disorder parameter T and doping. For the single crystal 
materials, one would identify the main origin of disorder with the change of Coulomb poten- 
tial at the Mn sites due to substitution of trivalent rare earths by divalent alkaline earths. 
This change amounts to some eV's, so that even including a reasonable screening it remains 
within the order of magnitude of the bandwidth. Crudely, we can assume that the different 
levels of doping changes the center of the distribution of energies but do not affect T. In 
that case, samples of Lai-^Sr^MnOs for example, differ only in the hole concentration x. 
According to our results shown in Fig. 3, the resistive behaviour of the samples with x=0.15, 
0.175, 0.2, 0.3 as reported in Ref. 16 by Tokura et. al. is consistent with a value of T of 
about 2 (in units of the hopping parameter t ~ 0.2eV). Cation or Oxigen vacancies should 
enhance the values of T favoring e localization. This is found in samples of (LaMn)i_ 7 3 
as reported in Ref. 17 and in samples of Lai.^Ca^MnO^, Ref 18. To be able to compare the 
temperature dependence or the resistivity with our results more presicely it is necessary to 
obtain the magnetization and the energy difference between the Fermi level and the mobility 
edge as functions of temperature. We will publish these results elsewhere. 

Other experimental results as transfer of optical weight ( that can be infered from Fig. 
2 ), or effect of pressure ( increase of t ), substitution of La by other rare earths ( decrease 
of t ), also find a simple and consistent explanation in terms of the model. As mentioned 
above, polaronic effects are expected to be present in these materials as in other perovskites 
in which electronic structure changes occour, but are not necessary to the understanding of 
the phenomena under consideration. 

In summary, we present here a simple model that can describe the properties of 
Lai_ x Aa;Mn03 and related compounds. It allows one to describe transitions between "metal- 
lic" and "insulating" states by shifts of the Fermi level relative to the mobility edge of the 
system as the magnetization changes. 
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FIGURE CAPTIONS 

Figure 1. Partial densities of states for spin up and down according to Eq. (6) for k — 5, 
E± = =f6, = 0, t — 1, and different values of the magnetization. 

Figure 2. Partial densities of states for spin up and down for the lower band including 
a Lorentzian distribution of diagonal energies of width r = 1, k — 5, t — 1, n — 0.88, and 
J » t. Bold lines indicate the zone of energies where the states are extended. The Fermi 
level (ep) is indicated for each value of the magnetization. 

Figure 3. Electron concentration versus magnetization "metal-insulator" phase diagram 
for k = 5, t = 1 and different values of T. 

Figure 4. Energy difference between paramagnetic and ferromagnetic phases in units of 



t for various values of T. 



